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Abstract. We derive a semiclassical theory for the detection of matter-waves. 
This theory draws on the theories of semiclassical optical detection and of fluid 
mechanics. We observe that the intrinsically dispersive nature of matter- waves is 
important in deriving such a theory. 



In nature particles exist in either of two forms: those with integer spin and those 
with half-integer spin. The former of these are known as bosons, and the latter as 
fcrmions. At a microscopic level, fermions are guided by the Pauli exclusion principle, 
which states that no two fermions can occupy the same quantum state. On the other 
hand, there is no such principle affecting bosons, and thus any number can occupy the 
same state. This is the basis of Bose-Einstein condensation, which was first proposed 
theoretically by Bose and Einstein in 1924 ^ E]. A Bose-Einstein condensate 
(BEC) is a system in which a macroscopic number of bosons occupy a single quantum 
state. To achieve this with a dilute gas of atomic bosons requires extremely low 
temperatures, such that the de-Broglie wavelength of the particles become larger than 
their mean spacing. Hence no BEC was experimentally realized in these systems until 
1995 ^flElEI using bosonic isotopes of Rb, Na and Li. 

A BEC is created by cooling the atoms using optical and magnetic forces, and 
then cooling them again using one of a number of techniques. The BEC is then held 
in a magnetic trap, which is switched off after a period of time to allow the atoms 
to expand so that imaging can take place. A recent experiment by Robert et al [7] 
involved the creation of a BEC of metastable triplet He (He*) and highlighted the 
ability to count single atoms falling from the trap after it was turned off. This allows 
for the exciting possibility of more detailed investigation of the quantum statistical 
properties of matter-waves. 

As a first step into this field, we must model the detection of matter-waves falling 
under gravity. We will use a simplistic model of a BEC, not including the effects of 
interactions which exist between atoms |H1 El • This will allow the features specific 
to detection to be more readily illustrated. As we are considering matter-waves, we 
can draw a direct analogy with the well known theory of the detection of light waves, 
or photo-detection, which will be outlined here. A more complete description can be 
found in many texts; e.g. |10|. 

The classical theory of photo-detection is based on the assumption that the 
probability of an ionization event occuring in the photo-detector in a time period 
dt is proportional to the cycle-averaged intensity I(t) of the incoming light: 

p(t)dt = £/(t)di, (1) 
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where £ is a constant of proportionality which represents the efficiency of the detector, 
including geometric factors such as its area, and dt is sufficiently small that the 
probability of more than one detection event occuring is negligible. In general, the 
cycle-averaged intensity is taken to be 

I(t) = le c\E(t)\ 2 ^cW(t), (2) 

where W(t) is the energy density. Under assumption if we take a time interval 
from t to t + T, then the probability of m detection events occuring is 

P m (T) = ( — exp[-fi]j, (3) 

where 

rt+T 

n = Zj dt'J(t') (4) 

and the angled brackets indicate a statistical average. From (j3J we can evaluate the 
mean number of detection events to be 

(m) = (n) . (5) 

We now wish to construct a semiclassical theory of matter-wave detection by 
analogy with the theory of photo-detection presented above. A natural way to proceed 
is to replace the electric field E(r, t) with the particle wavefunction ip(r, t). Thus the 
matter- wave analogy to the expression for I(t) in @ will be \ip(r, t)\ v, where we 
have included a characteristic velocity v. This is in direct analogy with the velocity 
of light c in the photo-detection theory and is of vectorial nature to allow for matter- 
waves which are not travelling perpendicular to the detector. It is also required so 
that the equations have the correct dimensionality. In the analysis that follows, v will 
be associated with the mean velocity of the wavepacket. The probability of detection 
over a time interval from t to t + T would again be given by and the average 
number of counts (m) by (0 where instead of (QJ , we have 

rt + T r 

n = Z J dt' J \ip(r,t')\ 2 v -dA. (6) 

We have now explicitly included the area of the detector A, and dA is the infinitesimal 
area element normal to the surface of the detector. If we assume that the particle 
wavefunction is normalized so that it contains on average N particles, then for all 
times t 

poo 

d 3 r \ip(r, t)\ 2 = N. (7) 

If the detector is of perfect efficiency then we would expect that for a wavepacket 
falling under gravity, a sufficiently long detection window and large detection area 
would produce a mean of N detection events. This means that from @ we might 
expect that as T — > oo, 

/OO P 
dt / mr,t)\ 2 v-dA = N, (8) 
-oo J A 

for the value of £ corresponding to a perfectly efficient detector. 

By drawing analogy with photo-detection of light waves, we have derived © 
which includes the characteristic velocity i). As a first approximation we might expect 
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that this will be the mean velocity of the wavepacket. This is not an approximation for 
light in free space because free space is not dispersive; at all frequencies light travels 
at c. For matter-waves, however, free space is dispersive. 

In order to take into account matter-wave dispersion we ought to base our theory 
of detection on the flux-density of particles: the mean rate at which particles cross a 
unit area of the detector. As particle number is a conserved quantity it must satisfy 
an equation of continuity [111 H^j 

J^(r,i)| 2 + V- J(r,t) = 0, (9) 

where J is the particle flux-density. This equation is of the same form as the one for 
local charge conservation in electromagnetic theory or, more relevantly for our purpose, 
relating particle density p and particle flux-density J = pv in fluid mechanics |13| . 
From (O we obtain a particle flux-density of the form 

J(r,t) = — Im {V*(r,t)W(r,t)}. (10) 

TO 

As this is analogous to the particle flux-density J = pv from fluid mechanics, it seems 
reasonable that © would become 

n = £^ At' J J(r,t')-dA. (11) 

Indeed, in electromagnetic theory a relation similar to © exists between the energy 
density Wit) and the Poynting vector (which gives the energy flux-density) |14l I15| . 
Equation is thus an approximation which holds in most experimentally realizable 
situations. In situations where this approximation is invalid, the cycle-averaged 
intensity in must be replaced by the magnitude of the Poynting vector. 

In order to illustrate fully the difference between the theories given by JfjJ) and 
(| 1 1 P - it is instructive to evaluate both expressions in the case of the detection of a 
wavepacket falling in the z-direction under gravity onto a flat, large- area detector 
aligned parallel to the x-y plane. Such a system closely models the He* experiment 
mentioned at the beginning of this paper, and it is one in which we would expect all 
particles to fall onto the detector, which will allow us to check the expression for (m). 

In evaluating the probability of detection for a wavepacket falling under gravity 
we will need to calculate the form of the matter- wave. We consider a model BEC, 
released at time t = 0, described by a Gaussian wavefunction centred at Tq with width 
parameter w, 

Mr,0) = A^(™ 2 )-3exp(-^— ^L\. (12) 

L 2w A J 

The standard solution to the Schrodinger equation takes the form 

^(r,*)=expJ^^U(r,0), (13) 



where H is the Hamiltonian, which in this case has the standard kinetic energy term 
and a gravitational potential term 
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We have taken the zero of gravitational potential energy to be at z = 0. When written 
in the position representation, this becomes 

H=-—V 2 + mgz, 15 
2m 

and so the particle wavefunction at a later time t will be 

i>(r,t) = A^5( 7ru; 2 ) -| exp |„| ^^-V 2 +mffz^|exp|- |r 2 J 2 o1 j. (16) 

Techniques outlined in ^B] allow us to evaluate the deriviatives in this expression and 
obtain a wavefunction for the matter-wave falling under gravity 

where we have defined the average "classical" position of the particle R(t) = (r(t)) = 
ro — ^gt 2 k, which gives the position of the centre of the wavepacket. As the wavepacket 
is accelerating from rest under gravity, the integral in JBJ will be given by 

J^(r,t)\ 2 v-dA = gt J J \rP(r,t)f dxdy. (18) 

It can be seen that the expression in l(T%|) depends on exp{ — (z— zq) 2 } and thus depends 
on the height that the wavepacket starts above the detection screen. With this taken 
into account, one can see that the integral of l|18[l over all time cannot give a constant 
value of N, and so the expression in (JSJ) cannot hold for any £ which is solely dependent 
on detector properties. This result can be verified numerically. 

If we now use (llOfl to calculate the flux-density of particles for this system, we 
obtain an expression for the integral in (|11|) 

J J \^(r,t)\ 2 dxdy. (19) 

It is straightforward to show that the integral of this expression over all time gives 
the average number of particles in the wavepacket N. Thus from (|llfl we can see that 
the constant of proportionality £ is in fact the efficiency of the detector rj, which takes 
values between and 1. 

It is clear to see that the expression obtained in (|19|) is that from l|18fl phis 
an additional correction, which is a height-dependent velocity term. This additional 
velocity term is a direct consequence of the dispersive nature of free space for matter- 
waves. From Ijl7|l it is clear that the wave undergoes dispersion as it falls under gravity. 
The detection theory based on (|18|l assumes that this dispersed wavepacket propagates 
through the detection plane at the mean packet velocity. The detection formula in 
<|19fl based on particle flux does not make this assumption and the factor %/m which 
quantifies the dispersion of the wave in H17[) also appears in the detection formula. If 
this factor is taken to zero either by taking K — > or m — ► oo, then the dispersion in 
<|17fl disappears, as does the additional velocity term in 119f) . The time variation of 
the integrals given by the two different theories are plotted in figure ^ where it can 
be seen that the differences in the expressions are quite pronounced: in a detection 
theory which takes account of dispersion the majority of particles will arrive earlier 
than they would in a detection theory in which dispersion is not correctly accounted 
for. 



/ J(r,t) 

J A 



■dA = 



gt- 



z-z a +gt 2 /2 
t + w 4 m 2 /(h 2 t) 
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Figure 1. Comparison of integrals in the intensity term for the correct and 
incorrect detection theories, as a function of time. The units have been chosen in 
such a way that g = h/m = 1, and we have chosen w = 1 and zo = 3. The solid 
line shows the correct theory 1191 and the dotted line shows the incorrect theory 

From l|17|) we can see that if the factor ht/(mw 2 ) is greater than unity, the 
wavepacket becomes significantly wider (due to dispersion) , and this dispersion ought 
to be taken account of in detection. As an example of how important dispersion is in 
the system under consideration, we take values from the He* experiment presented in 
. The time of flight of atoms here is 0.1s and the mass of a He* atom is 6.68 x 10 -27 
kg. We thus find that the dispersion factor will be important for any wavepacket with 
an initial width of less than 0.1mm. 

We have described in this paper the construction of a semiclassical theory of 
matter-wave detection, drawing on the well known theory of photo-detection. It is 
the intrinsically dispersive nature of matter-waves which prevents the direct analogy 
from working. We must instead consider the flux-density of particles, which gives an 
additional velocity term. Indeed if light passes through and is detected in a dispersive 
medium, the magnitude of the Poynting vector, which represents the flux-density of 
energy, must be used in place of ©• 

An instructive "next step" will be to consider the second quantized version of this 
theory. We would expect that in doing this, a situation of no detection in the early part 
of the wavepacket would feed back to modify the later part of the wavepacket. It is 
also clear that quantities other than particle number - such as energy, momentum and 
angular momentum - can be conserved. We intend to investigate these conservation 
laws, fluxes and the deposition of such quantities on a detection screen. 
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